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We present a prescription for generating pure spin current or spin selective current, based on 
adiabatic quantum pumping in a tight-binding model of a one dimensional conductor. A formula 
for the instantaneous pumped current is derived without introducing the scattering matrix. Our 
calculations indicate that some pumping cycles produce the maximum value 2 of pumped spin while 
others reverse the direction of current as a result of small alterations of the pumping cycle. We find 
pumping cycles which produce essentially any ratio of spin current to charge current. 
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I. INTRODUCTION 

The field of spintronics offers a vision of electronics 
that utilizes carrier spin in addition to carrier charge 1 . 
The rich potential of carrier spin for applications ranges 
from non-volatile devices to quantum computation. In 
order to realize this potential, however, it is essential to 
develop effective tools for the manipulation and trans- 
port of spin. Adiabatic quantum pumping is a mecha- 
nism that can transport charge at zero bias^ii^k. As a 
result of cycling two or more physical parameters that 
characterize a one-dimensional conductor, charge carri- 
ers get "swallowed" down the conductor like food down 
a throat, comprising a direct current. This method can 
deliver precise currents and requires no voltage bias. Re- 
cently, it has been shown that adiabatic quantum pump- 
ing in the presence of a magnetic field can also generate 
a spin currenii£i£i2i2iiS. For some fortuitous choices of ex- 
perimental parameters, it has even been possible to gen- 
erate a spin current without any charge current, which is 
termed a "pure" spin currenti£i2iifi. This is a promising 
development, and one wonders if it is possible to estab- 
lish complete control over both the amount of spin and 
the amount of charge pumped per cycle. 

A device for generating a pure spin current with im- 
proved control appears in Ref. |3, in which Zeeman energy 
is chosen as one of the adiabatic pumping parameters. 
There is no need to make a fortuitous choice of param- 
eters in this device - when the minimum and maximum 
Zeeman energies involved in the pumping cycle are equal 
in magnitude but opposite in sign, a pure spin current 
arises. However, if the maximum value of Zeeman en- 
ergy is not equal and opposite the minimum value, vari- 
ous combinations of spin current and charge current arise 
in a way that is not easily controlled. 

In this paper, we introduce a flexible approach to 
adiabatic pumping in which essentially any composition 
of spin/charge current can be chosen as desired. The 
scheme utilizes "generalized" pumping parameters each 
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of which depends on more than one physical parame- 
ter. With generalized pumping parameters, many differ- 
ent physical processes map to the same path in pumping 
parameter space. The result is greatly improved control 
over the pumping products. 

For instance, a pure spin current can be generated in 
the following way. In adiabatic quantum pumping, car- 
riers are transported with each cycling of the pumping 
parameters. By reversing the direction of the pumping 
cycle trajectory, one reverses the direction of the current 
flow. With generalized parameters, it is possible to make 
the spin-up pumping parameters traverse the exact same 
trajectory as the spin-down pumping parameters, but in 
the opposite direction. The result is that spin-up carri- 
ers get pumped in one direction and spin-down carriers 
in the other direction, leading to zero net transport of 
charge but non-zero spin current. This technique can 
be applied to any desired cycle in parameter space, so 
that the amount of spin pumped per cycle can be set by 
judicious choice of the trajectory in parameter space. 

Generalized pumping parameters also enable selective 
spin pumping wherein the current consists only of spin- up 
carriers or only of spin-down carriers. One spin's param- 
eters traverse a degenerate cycle that pumps no charge, 
while the other spin's parameters traverse a productive 
cycle. This selective spin pumping is a valuable tool; 
by combining and repeating spin selective pumping, it 
is possible to generate any rational proportion of spin 
current to charge current. 

Finally, in addition to pure spin pumping and spin se- 
lective pumping, we consider a family of cycles which 
produce arbitrary ratios of spin current to charge current 
after exactly one cycle. Unlike the schemes mentioned 
above, the correct cycle in this case cannot be fixed in a 
deterministic way; trial and error is necessary. However, 
we present an argument that an appropriate cycle gener- 
ally exists. Moreover given a cycle that produces a given 
ratio of spin current to charge current, we also show how 
to traverse a cycle that produces the inverse ratio. 

The rest of the paper includes a presentation of our 
model Hamiltonian, in sec. II, and a derivation of the 
pumped current. Section III. discusses the generation of 
a pure spin current. In Sec. IV, we present a method 
for selective spin pumping of just one spin orientation. A 
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means of generating arbitrary ratios of spin current and 
charge current is proposed in Sec. V, and we conclude 
with Sec. VI. 



II. MODEL HAMILTONIAN AND PUMPED 
CURRENT 



We consider transport through a one-dimensional 
channel of sites, schematically depicted in Fig. [I] The 
following is our model Hamiltonian 

H = + + 
Ho = -JJ2 

a \i+la a n<J + a n(j a n+lcr 

Vi = ^ u-i n-icj + m nig 

a 

V 2 = ^-aE^n-ir-aE? n la . (1) 



In the Hamiltonian,— J is the nearest neighbor hopping 
amplitude, o) no . is the electron creation operator at site 
n for spin a and ni a = aj a aia is the number operator 
at site I for spin a. The first term Ho is the Hamilto- 
nian of carriers in a homogeneous channel and V\ is the 
impurity potential. The impurities are characterized by 
on-site energies {u-i,ui}, which simulate potential bar- 
riers at sites —I and I. The V2 term contains the Zeeman 
spin energies E^_ x = gusH-h Ef = gfisBi for carriers 
in the localized external magnetic fields at sites —I and I. 
We assume that all four experimental parameters in Fig. 
[I] can be tuned precisely; naturally, this would be chal- 
lenging to realize in the laboratory. We set spin |a| equal 
to 1 instead of 1/2. We assume that the four quantities 
{u-i, ui, E^_ x , E^} are slowly varying time-dependent pa- 
rameters. We define generalized spin-dependent pumping 
parameters Xi a }. 
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FIG. 1: Gray dots located at —I and / sites represent impu- 
rities on a wire; the impurities are characterized by potential 
barriers u-i,ui. Localized magnetic fields B-i,Bi are applied 
on those impurities. 



in terms of which the Hamiltonian becomes 
H = H + V 

Ho = -J ^ cJ n+la a nc7 + Q^Qn+lo 

71,(7 

V = J ^ X ~icr n ~lcr + UUj. 



(3) 



For mathematical convenience the following are set to 
unity Ti = e = lattice spacing = 1. By defining the spin- 
dependent parameters we make explicit the 
fact that spin-up carriers and spin-down carriers are sep- 
arately controllable. 

Instantaneous scattering states \Xpa) of the Hamilto- 
nian are obtained by ignoring the time-dependence of the 
pumping parameters and using the Lippmann-Schwinger 
equation 
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\ Xp a) = [l + G(E p )V}cl <7 \0) 



(4) 



where cj a is the carrier creation operator for energy 
E p = — 2 J cos p and spina, and G(E p ) = l/{E p — H+irj) 
is the retarded full Green's function, 77 = + . The time- 
variation of the potential is adiabatic if it is slow com- 
pared to the dwell time of a carrier in the scattering 
region^. The time-dependence is then restored to first 
order by adiabatic corrections 4 



IrV) = \x P a) -iG(E p ) \x pa ) 



(5) 



In terms of these first order scattering states JJJ, the 
instantaneous pumped current associated with spin a is 



a(n) = I" dE F(E) \-2J J S(E - E p ) x 



(na\(/) pa )(n + lcr^)* 



(6) 



where ^5 indicates imaginary part, F(E) is the Fermi dis- 
tribution function, and \no) = a^O). 
First the integral over p is evaluated 



2J3 ^ ^6(E-E p )(n*\<f> pa )(n + l<T\(l> P a)* (7) 
2 J 2 9 V X ma r 6(E - E p ) x 

i{n + la\ G 2 (E P ) |raa)*(na|xpa)(rncr|xp<r)* 



-i(na\ G 2 (E P ) \ma) {ma\x p( r) (n + 1(j\x P ctY 



(8) 



where the identity \x p(J ) = G(E p )V\x P a) has been used 4 . 
Eq. JSJ is evaluated using the identity 
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G(ncr, ma; E) — G* (ncr, ma; E) 



(9) 
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where G(na, racr; E) = (na\G(E)\ma). The result is 



where 



(n + lcr\G 2 (E)\ma)* x 



j2 

^5 / Xracr 
m—±l 

[G(na, ma; E) — G* (na, ma; E)] 
— {na \G 2 (E)\ma) x 

[G(ma, n + la; E) — G* (ma, n + la; E) 



(10) 



Eq. can be simplified if we use the fact that the 
1-D Green's function takes a plane- wave form at large 
distances^ 



G(n + la, ma; Ek) 



e tk G(na,ma; Ek) for n — > oo. 

(ii) 



The asymptotic condition n — > oo means the observation 
point is far away from the scattering center and is located 
on the right side. Inserting this into Eq. ([TUj) produces 
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— ^ Xma dE[G(na,ma] E)G* (n + la,ma; E)]. 

m—^zl 

(12) 

At the last step of the calculation, we use the identity 
(na\G 2 (E)\ma) = — <9#G(na, raa; £?). This identity is 
derived by expanding G(E) in a basis of energy eigen- 

states of H, G(E) = ^ ^-e^ ' The integration over 
energies in Eq. (0) can be performed with the result ([T2j) 
at zero temperature where the Fermi distribution func- 
tion is a step function 



J 2 

/ Xma 

m=±l 



^s\G(na, ma] E)G* (n + la, ma; E)~\ \e=e f for n — > oo. 

(13) 

This closed form for the instantaneous pumped current 
in terms of retarded full Green's functions is one of the 
main results of this paper. In contrast to most treatments 
of quantum pumping, our derivation does not introduce 
the scattering matrix. (However, see e.g. Ref. 13.) The 
identity © is crucial in our derivation. 

To evaluate the pumped current associated with the 
model Hamiltonian JJJ, we need to compute the full 
Green's function. The full Green's function can be ex- 
pressed in terms of the free Green's function by the al- 
gebraic identity 11 G(E) = G (E) + G(E)VG (E), where 
G (E) = 1/(E - Ho + irj). Eq. O is evaluated to be 



ja(n) 



= ±l 



; [f(-ma)f(-ma) + h(-ma)h*(-ma)}%[G o (0a, 0a; E)] 



+M[f(-ma)h*(-ma)]%[Go(ma, -ma;E)] 
—sign(m)^s[f(—ma)h*(—ma)]R[Go(ma, —ma; E)] 



(14) 
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J[X- ma + X m *]G o (0*, 0a; E) 



J cJ 2 X-maXmaG 2 ) (0a, 0a; E) 



-J 2 X- 



7 X ma G (ma, —ma; E). 



Eq. (T2J) can be evaluated since the explicit expression 
for the free Green's function 11 is just Go(na^ma; Ek) = 
e ik\n-m\ Js'mk. The pumped charge associated with 
spin o~ after one pumping cycle is obtained by integrating 
Eq. fnj with respect to time 



(15) 



q a = (b dtj a (ri). 



The pumped charge associated with spin a (|T5|) can 
be represented as a surface integral rather than a line 
integral 2 



Qa — — / [ dX-iadXia sign{m) x 

71 J Js jrl, 

d Xrn(T ^[G(na, 



m=±Z 

ma; E)G*(n + la, —ma; E)] 1^=^ 



(16) 



where S indicates the area which is enclosed by the 
pumping cycle in parameter space. There are two fea- 
tures of the pumped entity (charge or spin), evident in 
the last two equations, that we utilize to control the flow 
of spin and charge. First, in the line integral form in 
Eq. (|T5J) it is clear that reversing the direction of the time 
cycle changes the sign of the integral, indicating a flow in 
the opposite direction over a pump cycle. Secondly, the 
surface integral form, Eq. ([TB|) . shows that the magnitude 
of the pumped quantity in a full cycle depends entirely 
on the enclosed surface S in parameter space. Therefore, 
surfaces that enclose areas with identical functional form 
and values of the integrand will yield identical magni- 
tudes of the pumped charge or spin, while the direction 
of traversal of the bounding curve will determine the di- 
rection of flow. 



III. SPIN CURRENT WITHOUT CHARGE 
CURRENT 

In the previous section we established how we can con- 
trol the direction and magnitude of the pumped charge 
associated with each spin state, and found expressions 
to determine them. In particular, varying the mag- 
netic fields in the generalized parameters allows differ- 
ential manipulations of up and down spin states, because 
the path in parameter space of the spin-up parameters 
{X_/|,X/|} becomes distinct from the path in parame- 
ter space of the spin-down parameters We 
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will now use these considerations to present two dis- 
tinct types of pumping cycles which generate only a pure 
spin current, with zero transported charge after each cy- 
cle of pumping. The first cycle relies on the fact that 
the integrand of ((16)1 is symmetric under exchange of 
X-i a and Xi a , as one expects given the form of the 
Hamiltonian (0. Consider two square cycles (we will 
use the term "boxes") of side length 5 which are lo- 
cated symmetrically in parameter space about the line 
X-i — X\ (see FigEt)- For the cycle taken by spin- 
down parameters {X-i^X^}, we pick an arbitrary point 
(a, b)^ as the initial choice of parameters. For the cy- 
cle taken by spin-up parameters {X-i^X^}, the ini- 
tial point (6, a — 5)^ is chosen. For those initial points, 
we find U-i/J = (a + 6)/2, u\j J = (a + b — 5)/2 and 
E z _ x /J = (a - 6)/2, Ef '/J =-(a-b- 5)/2. We fix u- h 
Ef throughout the pumping cycle. First, the Zeeman 
energy at site —I divided by J, E Z J J, is decreased by 5, 
from (a — b)/2 to (a — b — 25)/ 2. The resulting motion 
in parameter space is parallel to the X-\ axis, with the 
spin-up parameters moving in the positive direction and 
the spin-down parameters moving in the negative direc- 
tion. Next, the potential barrier at site I divided by J, 
u\j J, is increased by the amount of 5, from (a + b — S)/2 
to (a + b + S)/2. The spin-down and spin- up parameters 
both shift upward parallel to the X\ axis. Next, the Zee- 
man energy E^_J J is increased by 5, from (a — b — 28)/ 2 
back to {a — b)/ 2. Finally, the potential barrier u\j J is 
decreased by 5, from (a + b + 5)/2 back to (a-\-b — S)/2, to 
complete the cycle. The form of the definition ensures 
that spin-up and spin-down parameters shift in opposite 
directions when the Zeeman energy is varied and shift 
in the same direction when the potential barrier, which 
results from an electrical potential, is varied. The combi- 
nation of these two effects moves the spin-up parameters 
in a counterclockwise cycle, and the spin-down parame- 
ters in a clockwise cycle. In addition, our chosen steps 
generate two square cycles located symmetrically about 
the line X-\ = X\ in parameter space. Because of the 
symmetry in (JT5)k these cycles lead to zero total pumped 
charge q c = + q± = 0. On the other hand, the pumped 
spin is q s = — q± ^ as long as q^(= q\) ^ 0. The 
result is a pure spin current. 

A second type of cycle generates a pure spin cur- 
rent without relying on the symmetry between X\ and 
X-i. Consider a rectangular box in parameter space (see 
FigEb). By choosing two initial points appropriately, 
we can make the parameters execute cycles on the same 
rectangular box of width 5 and height 7 but in opposite 
directions. For the cycle of spin- up parameters, we pick 
an arbitrary point in parameter space (a, For the cy- 
cle of spin-down parameters, the initial point (a + 5, b)i 
is chosen. These choices correspond to U-i/J = a + 5/2, 
ui/J = 6, E Z JJ = 5/2, and Ef /J = 0. We vary u t 
and E Z L while fixing U-i and Ef . First the Zeeman en- 
ergy E^J J is decreased by the amount 5, from 8/2 to 
—6/2. Second, the potential barrier ui/J is increased by 
7, from b to b + 7. Third, the Zeeman energy E^J J is 




FIG. 2: (a) Two identical square boxes are located symmet- 
rically about the X-i = Xi line. For spin up, the parameters 
{X-i^Xi^} begin at (6, a — S) and go counterclockwise. For 
spin down, the {X-ii,Xii} begin at (a, b) and go clockwise, 
(b) The spin up and spin down cycles are located on the same 
rectangular box. For spin up, the cycle begins at (a, b) and 
goes counterclockwise. For spin down, the cycle begins at 
(a + 8, b) and goes clockwise. 



increased by 5, from —5/2 back to 5/2. Finally, the po- 
tential barrier u\/ J is decreased by 7, from 6 + 7 back to 
b. As a result of these variations, the spin-up parameters 
traverse the rectangular box in the counterclockwise di- 
rection while the spin-down parameters traverse the rect- 
angular box in the clockwise direction. Since the cycles 
enclose the same region, but move in opposite directions, 
a pure spin current arises. 

In this second type of cycle, note that Ef is fixed at 
zero throughout the pumping. This suggests a means of 
realizing the cycle experimentally. Rather than trying to 
produce a localized magnetic field B_i, one could apply 
a global magnetic field. If all sites except for the site —I 
have a negligible g-factor, the desired Hamiltonian Q 
will arise—. 

We conclude based on the above analysis that, for 
any given pair of identical symmetrically located square 
boxes or for any given single rectangular box in param- 
eter space, there always exists a pumping cycle which 
generates a pure spin current. This finding implies great 
flexibility in the control of pumped pure spin after one 
cycle. Since the quantity of pumped spin depends on the 
shape of the enclosed area and its location in parame- 
ter space, one can tune the quantity of pumped spin by 
changing these characteristics of the pumping cycle. 

The following plots made using the expressions derived 
in Sec. II demonstrate flexibility in controlling a pure 
spin current. Fig. 03 and Fig. EI portray how the total 
pumped spin depends on the size of the box in parame- 
ter space and on its location. For the first type of cycle 
involving pairs of symmetrically positioned square boxes 
in parameter space, three different sizes for the box pairs 
are considered in Fig. [HK- All meet at the same point on 
the line X_i = X\. The dependence of the pumped spin 
on the location of the cycle in parameter space is studied 
by moving that common meeting point a distance \[2d 
along the line X_\ = X\. Each curve in Fig. shows 
the variation of the pumped spin as a function of d for a 
specific box size, a fixed Fermi wave vector kp = 1.4 and 
impurities at ±1 = ±1. The different curves correspond 
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FIG. 3: (a) Three pairs of square boxes are located symmetri- 
cally about the median X-i = Xi and share the same meeting 
point. We move those pairs along the X-i = Xi line. For the 
three pairs, the box side lengths are 2,4, and 6 respectively. 
The pairs shift a distance \[2d along the X-i = Xi line, (b) 
Pumped spin q s vs d at fixed k F — 1.4 and for 1 = 1: The solid 
line is for the pair of large boxes of side length 6, the long- 
dashed line is for the pair of medium boxes of side length 4, 
and the short-dashed line is for the pair of small boxes of side 
length 2. The result shows monotonic increase of pumped 
spin with box size and monotonic decrease of pumped spin 
with projected distance d. 



to the three different box sizes. The plots show a mono- 
tonic increase of pumped spin in a cycle as the box size 
increases (except where the pumped spin vanishes for all 
three box sizes). This is physically reasonable since the 
box size determines the difference between the minimum 
and maximum value of the potential barriers. For large 
boxes, the potential barriers change a lot during the cy- 
cle, resulting in more pumped current, and the opposite 
is true for small boxes. The pumped spin decreases as 
the distance d increases. As d increases, the minimum 
height attained by the potential barriers gets larger. As 
a result, the current must traverse an increased potential 
barrier, so that the transmission is decreased. 

For the second type of cycle discussed above for gener- 
ating a pure spin current, we consider three square boxes 
that each straddle the median line X-i = X\ as shown in 
Fig. The pumped spin in a cycle is plotted as a func- 
tion of d where the lower left corner of the box is at the 
point (d, d). The two plots, Fig. and Fig. corre- 
spond to different choices of Fermi wave vector, kp = 1.4 
and kp = 3.1 respectively; the clear differences between 
the two sets of curves indicate that the Fermi wave vector 
is another essential factor in controlling the spin current 
flow. For both plots the impurity locations are ±Z = ±1. 
In Fig. ^3, the pumped spin increases monotonically with 
box size at most locations d. For large d, the behav- 
ior of the curves in Fig. 2t> shows monotonic decrease 
with d like the curves in Fig. [SJd. However, for small d 
the pumped spin increases with d, so that each curve in 
Fig. 0Jd has maximum pumped spin between d = 3 and 
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FIG. 4: Plots of q s vs d. In each plot, the solid line is for 
a large box of side length 6, the long-dashed line is for a 
medium box of side length 4, and the short-dashed line is for 
a small box of side length 2. (a) We move each square box 
along the line X-i = Xi, where the position of the lower left 
corner is (d,d). (b) Plot at k F = 1.4 and for / = 1. For small 
d, negative pumped spin is found for the smaller two boxes 
(i.e. spin is pumped in the opposite direction), (c) Plot at 
k F = 3.1 and for / = 1. The maximum value 2 of pumped 
spin per cycle occurs on the solid line. The pumped spin q s 
is independent of box size for sufficiently large d. 



d = 5. This kind of maximum can be explained in terms 
of resonant transmission 15 . When the pumping cycle in- 
cludes locations in parameter space for which the Fermi 
energy satisfies a resonance condition, an enhanced trans- 
mission coefficient leads to a large pumped current. For 
the parameters chosen in the figure, a line of resonant 
points in parameter space runs near (X^X-i) = (6,6). 
As each box shifts with increasing d to enclose these 
resonant points, the pumped current grows even though 
the minimum heights of the potential barriers get larger. 
For d near zero in the case of the two smaller boxes in 
Fig. ^3, there is actually a region of negative values for 
the pumped spin, indicating a reversal of direction of the 
spin flow. The direction of pumped spin changes as the 
box moves or its size increases. This shows we can con- 
trol the direction of the spin current without reversing 
the whole cycle, by simply adjusting the box size or lo- 
cation. 

For Fig. the maximum pumped spin for the largest 
box is 2, so that during each cycle one spin- up carrier 
goes to the right and one spin-down carrier goes to the 
left. This large maximum value is due to both the low 
minimum barrier heights and the resonance transmission. 
For kp = 3.1, a resonant line runs near (X^X-i) = (1,1) 
and boxes enclosing the resonant line also have the low 
minimum barrier heights, by adding up two effects, the 
pumped spin has the large maximum value such as 2. 
With increasing d from the origin, the pumped spin de- 
creases to assume the same finite value for the three dif- 
ferent box sizes, indicating that there is little variation in 
the integrand in Eq. (|TfiJ) in the surface integral far from 
the origin whereas there are stronger variations close to 
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FIG. 5: (a) An example of pure spin current. There is no net 
flow of charge, (b) An example of selective spin current. Only 
one kind of spin contributes to current, there is no cancellation 
of charge. 
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FIG. 6: (a) Pumping cycle for spin-up parameters, beginning 
at the point (a,b) in the generalized parameter space and 
proceeding counterclockwise. It encloses a non-zero area, and 
generates a non-zero spin- up current, (b) Pumping cycle for 
spin-down parameters, also beginning at (a,b). This cycle 
encloses no area, so there is no spin-down current generated. 



the origin. 



IV. SELECTIVE SPIN PUMPING 

We consider another type of cycle which has a spin- 
filtering effect. This cycle selectively pumps one kind of 
spin, so that an equal spin and charge current flow (See 
Fig. 0). Consider a rectangular box in parameter space. 
We choose the same initial point (a, b)^i for both spin- up 
and spin-down parameters. This choice of initial point 
implies that U-i/J = a,ui/J = b and E^_ t = Ef = 0. 
We fix Ef at zero. First, we increase U-\ and decrease 
simultaneously in such a fashion that (u-i + E^_ t )/ J 
remains at the initial value a while (u-i — E^)/ J is in- 
creased by the amount c), from a to a + 5. At the end of 
this process, U-i/J is a + 8/2 and E^J J is —5/2. The 
spin-up parameters shift parallel to but the spin- 

down parameters remain unchanged. Second, u\j J is in- 
creased by 7, from b to b + 7. The spin-down parameters 
and the spin-up parameters both shift parallel to the X\ 
axis. Next, we decrease u-\ and increase E^_ x simultane- 
ously, keeping U-i+E^_ l fixed while (ui—Ef)/ J decreases 
by <S, from a + S back to a. This produces a path paral- 
lel to X-i for spin- up parameters, but does not shift the 
spin-down parameters at all. Finally, u\/ J is decreased 
by 7, from 6+7 back to 6, shifting both spin-down param- 
eters and spin- up parameters parallel to X\ . For spin- up 
parameters, this cycle makes rectangular path enclosing 
a non-zero area (see Fig. that produces a current of 
spin-up carriers. On the other hand, the cycle for spin- 
down shifts along a straight line in parameter space that 
encloses no area (see Fig. EJJd) and pumps no current. 
The result is perfect selective spin pumping of spin-up 
carriers. Naturally, a selective current of spin-down car- 
riers can be generated with trivial modifications to this 
protocol. 

More generally, we can transfer charge and spin to 
achieve any rational value of q s /q c by combining and 
repeating spin selective cycles. Suppose that the value 
Qs/Qc = M/N is desired, where M and N are integers. It 



is always possible to find two integers n and m satisfying 
n/m = (N — M)/(N + M). By performing \m\ selective 
spin- up cycles and \n\ selective spin-down cycles, we can 
generate an arbitrary rational value for q s /q c - (If m is 
positive, the spin-up cycles should be traversed in a coun- 
terclockwise direction, while for negative m they should 
be traversed in a clockwise direction. The same is true 
for n and the spin-down cycles.) 



V. ARBITRARY COMBINATIONS OF SPIN 
CURRENT AND CHARGE CURRENT 

In earlier sections, we gave definite cycles that could 
be used to pump spin with no charge, to selectively pump 
carriers of a given spin orientation, or to pump a rational 
ratio of spin to charge. Here, we argue that other cycles 
can produce arbitrary ratios of spin current to charge 
current, requiring only one cycle of pumping with no rep- 
etition. (However, we do not give a recipe for identifying 
the cycle; trial and error tuning may be needed.) 

Consider a single rectangular box and put two congru- 
ent rectangular cycles symmetrically at its ends as shown 
in Fig. EJl. The left cycle is traversed by the spin-up 
parameters, and the right cycle by the spin-down param- 
eters. We first describe a protocol in which the spin- 
up and spin-down parameter cycles are both traversed 
in a counterclockwise direction. Choose initial points 
at the lower left corner of each cycle. For spin-up pa- 
rameters, the point is (a, 6)j. For spin-down parameters, 
the initial point is (a + (3 — J, b)^. These choices imply 
that the physical parameters are u-i/ J = a + (J3 — 5)/2, 
ui/J = 6, and E^/J ={(3- S)/2, and Ef = 0. Suppose 
that the Zeeman energies E* x and Ef are fixed. As a 
first step, we increase U-i/J by 5, from a + (/? — S)/2 to 
a + ((3 + S)/2. Next, we increase u\/ J by 7, from b to 
6 + 7. Then, we decrease U-\j J by (5, from a + ((3 + S)/2 
back to a + ((3 — S)/2. Finally, we decrease u\j J by 7, 
from 6 + 7 back to b. This is the complete protocol. For 
counterclockwise traversal, it is typically the case that 
the charge produced by the spin-up cycle will be posi- 
tive and so will the charge q^ produced by the spin-down 
cycle. As a result, the ratio q s /q c = (#t — #|)/(#T wm 
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FIG. 7: (a) We consider two congruent rectangular boxes 
symmetrically located within one large rectangular box. By 
varying 8, we shift the interior edges of the rectangular boxes. 
The left rectangular box shows the trajectory of spin-up 
pumping parameters, while the right rectangular box is for 
spin-down pumping parameters. For the first type of cycle 
(shown here), both the small rectangles are traversed in a 
counterclockwise direction. For the second type of cycle (not 
shown here), spin-up parameters traverse in a counterclock- 
wise direction while spin-down parameters traverse in a clock- 
wise direction, (b) Plot of q s /q c vs 5 for first type of cycle at 
k F 3.1 and for I = 1. We set (a, b) (0, 0) and (3 = 7 = 8. 

5 is varied from 1 to 8. Any q s /q c in the range to 1 can be 
found. Note that the origin of plot is (1,0), not (0,0). 

typically satisfy \q 3 /q c \ < 1- Tuning the parameters will 
typically permit any desired value of the ratio, as shown 
in Fig. 03. In order to attain a ratio \q 3 /q c \ > 1, one can 
use the same two rectangular cycles, traversed in, say, a 
counterclockwise direction for spin-up and in a clockwise 
direction for spin-down. The lower left point in the spin- 
up rectangle and the lower right point in the spin-down 
rectangle serve as initial points. It follows that the phys- 
ical parameters take the initial values U-i/J — a + (3/2, 
ml J = b, E z _ x /J = and E? = 0. We fix u_ u E? 
during the cycle. First, E^JJ is decreased by 5, from 
(3/2 to (3/2 — 5. Next, u\/ J is increased by 7, from b to 

6 + 7. Third, E^J J is increased by 5, from (3/2 — 5 back 
to (3/2. Finally, u\/ J is decreased by 7, from 6 + 7 back 
to b. Since we are traversing the same two rectangles as 
the previous protocol, we see that the same value of 
will be produced, but the spin-down charge will now be 
—qi, where q^ is defined as the spin-down charge pro- 



duce by a counterclockwise traversal. The result is that 
Qs/Qc = (<?T + Ql)/(Q} ~ qi)i which is simply the inverse 
of the value of the ratio obtained in the first protocol, so 
that now \q s /q c \ > 1 typically. Given these two proto- 
cols, we should be able to achieve arbitrary combinations 
of spin current and charge current over the whole range 
< \q 3 /q c \ < +°° by varying 6 (see Fig [3d) which de- 
termines each pumping cycle within the rectangular box. 
Although tuning is required to achieve a given arbitrary 
ratio, the exact inverse ratio for that combination can 
be attained predictably by reversing one of cycles. The 
extreme cases \q 3 /q c \ = +00 and \q s /q c \ — 1 correspond 
respectively to the cases of pure spin pumping and to 
spin selective pumping described above. If q s /q c is posi- 
tive, a corresponding negative ratio, which has the same 
absolute value, can be obtained by exchanging cycles so 
that the left box is for the spin down parameters and the 
right box is for the spin up parameters. Thus, all possible 
ratios q s /q c are attainable. 



VI. CONCLUSION 

In conclusion, we have presented a deterministic way 
to produce a pure spin current, a spin selective current, 
or a rational ratio of spin to charge current. Our proposal 
relies on generalized pumping parameters, each of which 
depends on more than one physical parameter. In our 
calculations, the maximum value 2 for the pumped spin 
is observed for some pumping cycles, and we find that 
the direction of the spin current can be manipulated via 
the size or location of the pumping cycle in parameter 
space. We also presented an argument to show that it is 
typically possible in a single cycle to pump an arbitrary 
ratio of spin current to charge current (although some 
trial and error may be needed to find the right cycle). 
These results suggest that adiabatic quantum pumping 
could be a versatile tool for generating a desired current 
in a spintronics device. 
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